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Abstract: Several solenoids are usually installed in electron cooler device to guide the motion of the electron beam 
in the cooler. However, the solenoids also have influence to the ion beam in the cooler storage ring. The transverse 
motion of the ion beam in storage ring will become coupled, if the solenoids installed in the electron cooler are not 
compensated perfectly. In this paper, the coupled transverse motion due to the uncompensated cooler’s solenoids 
of CSRm (The main storage ring in the IMP, Lan Zhou, China) is studied, and the coupled beam’s envelopes are 
calculated by a new method. 
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1 Introduction 

An electron cooler device named EC-35 is installed 
in the CSRm to improve the quality of the ion beam in 
the storage ring[l]. As Figure 1 shows, there are several 
solenoids installed in the electron cooler, and revolution 
ion beam passes through the whole main solenoid and 
parts of the two toroid every turn. If the solenoids are 
not compensated, it will make the beam’s betatron mo¬ 
tion of horizontal and vertical directions become coupled. 
Usually, two compensation solenoids installed outside the 
electron cooler are used to deconole the beam. 



Fig. 1. The general arrangement of the EC-35, 1, 

3 are the toroids, 2 is the main solenoid. 

In this paper the transverse motion of ion beam in 
CSRm will be studied when the compensation solenoids 


are not used or are turned off suddenly. To research 
the coupled motion two methods are used: the tracking 
method and the Ripken’s analytical method. 

In Section 2 the simulation results of the coupled be¬ 
tatron motion will be presented, in Section 3 the Rip¬ 
ken’s method is used to parameterize two dimensional 
coupled betatron motion of CSRm, and the boundary 
functions in phase plane for coupled single particle and 
beam are obtained. In Section 4 a new method used 
to calculate the coupled beam envelopes are introduced 
and the calculating results are presented. In Section 5 
the conclusions are summarized. 


2 Simulation of the coupled betatron 
motion 


The 4*4 transmission matrixes of every element in 
the ring CSRm are used to simulate the coupled beta¬ 
tron motion. Particularly, transmission matrixes of the 
solenoid are presented by Equation (1), when the direc¬ 
tion of the longitudinal magnetic fields in the solenoid is 
opposite to the particles’ motion direction[2] (The situa¬ 
tion in the CSRm): 
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value of longitudinal fields of solenoid, L is the length of 
solenoid, G is the magnetic rigidity of ion particle. 



Fig. 2. The longitudinal component of magnet 
fields in the cooler section of CSRm 

In our case, the longitudinal magnetic fields in the 
cooler section of CSRm are not uniform, the value 
changes along longitudinal position, as Figure 2 shows. 
So we can not treat the whole cooler section as a solenoid. 
In order to obtain more accurate results, the transmis¬ 
sion matrixes of the solenoids in cooler section are ob¬ 
tained by following ways: 

Divide the cooler section into several parts, each part 
acting as a independent solenoid has a 4*4 transmission 
matrixes. By averaging the magnetic data of each part, 
the Bz can be obtained, then using equation (1) the 
transmission matrixes of the each part can be calculated 
out. 

7MeV/u 12 ( 76 + -g cPosen as example particle in simu¬ 
lation. Assuming that the initial particles matching the 
uncoupled twiss function at the starting point, and the 
emittances are equal to the acceptances of the CSRm. 
The detail parameters for initial particles are shown in 
Table 1[3]. 


Table 1. Parameters used in the calculations. 


Particle 

7MeV/u 12(76+ 

Horizontal acceptance Ax 

2007rmmmra,(i 

Vertical acceptance Ay 

^OTTmmmrad 

Horizontal {3 function 

I3x = 4.68m 

Vertical ^ function 

I3y = 30.17m 

Horizontal a function 

P 

II 

o 

Vertical a function 

P 

II 

o 


1 > Single particle’s motion 

A particle is generated randomly in (x, x’, y, y’) phase 
plane, then use transmission matrixes to track the par¬ 
ticle turn by turn. At one fixed watching point, the par¬ 
ticle’s projection in phase plane can be obtained turn by 
turn, as Figure 3,4 show. 



-2.51- ^^^^^ -' 

-0.03 -0.02 -0.01 0 0.01 0.02 0.03 

X [m] 


Fig. 3. The single particle’s projection in (x,x’) 
phase plane for 20000 turns in one watching point 
of the ring. Red point is the initial position of the 
particle. 

Form the Figure 3, one can find that: with coupling, 
a single particle’s projection for 20000 turns is different 
from the uncoupled case, the coupled particle’s projec¬ 
tion for 20000 turns at watching point no longer be a 
ellipse, instead it’s a torus. 



-81 -^^^-, 

-0.1 -0.05 0 0.05 0.1 

X [m] 


Fig. 4. Two particles’ projection in (x,x’) phase 
plane for 2000 turns in a watching point of the 
ring. The red point are for one particle, the blue 
points are for another. 

Further simulation indicates that: 1. For a same par¬ 
ticle, the shape of the projection torus is different at dif¬ 
ferent watching points of the ring. 2. For different initial 
particles, the shapes of the projection torus’s boundary 
are same at same watching poind, but the sizes are differ¬ 
ent (Figure 4 shows). 3. For a particle the boundary of 
the projection torus at every watching point will become 
fixed after several turns, the particle’s footprint will not 
beyond the boundary with the increase of the turn. 

2 > Beam’s motion 
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To simulate the beam’s motion in coupled situation, 
1000 particles are generated with the parameters shown 
in Table 1. Then use the transmission matrixes to track 
the particles turn by turn. The 1000 particles’ footprints 
in phase plane at one watching point for 2000 turns are 
shown in Figure 5. 



.gi-1-,-^^^^-, 

-0.2 -0.15 -0.1 -0.05 0 0.05 0.1 0.15 

y [m] 

Fig. 5. Beam’s projection in (y,y’) phase plane for 
2000 turns in a watching point of the ring. The 
yellow ellipse is the initial matched emittance cir¬ 
cle. Red points are the initial beam. 

From Figure 5 one can find that: with coupling, the 
beam’s projection for 2000 turns in phase plane is be¬ 
yond the initial emittance circle. Same as the uncoupled 
casethe boundary of the projection is also a ellipse. 

Further simulation indicates that: 1. The boundary 
of the beam’s projection in phase plane is different ellipse 
at different watching points. 2. For certain watching 
point the boundary ellipse is fixed, and not changes with 
the revolution of the beam. In other words, the envelopes 
along the ring for certain initial beam are certain, just 
like the uncoupled case. 

3 > The variation of emittance 

The variation of emittances in horizontal and in verti¬ 
cal directions are also calculated turn by turn, the result 
is shown in Figure 6. 



Fig. 6. The emittances’ changes turn by turn. 


From Figure 6, it can be found that the emittances 
in horizontal and vertical directions are not longer con¬ 
stants, however, the sum of the two emittances is a in¬ 
variant. 

From above simulation results, it become clear that 
the coupled betatron motion of beam is quite different 
from the uncoupled case. In coupled case, the beam 
emittances is no longer invariant. Although the beam 
envelope is existing, the equation r = y/e^ parameterized 
by Courant-Snyder parameters is no longer available to 
calculate the couple beam’s envelope. To calculate the 
coupled beam envelope and represent the coupled beam 
motion of CSRm, the method parameterize two dimen¬ 
sional linearly coupled betatron motion should be intro¬ 
duced. 

3 The coupled motion’s parameteriza- 
tions 

There are several methods usually used to parame¬ 
terize the coupled betatron motion. Ripken’s approach 
and the approach obtained by Edwards and Teng are 
most well-known of them[4]. Although the Ripken’s ap¬ 
proach is more complicated, it can give straightforward 
expressions of the particle’s coordinates in the laboratory 
frame[5]. So Ripken’s approach is chosen to parameter¬ 
ize our coupling situation. In Ripken’s theory a particle’s 
trajectory can be expressed by several lattice functions 
as [5]: 

/ X \ 

x' 

y 

\ y' J 

( cos (4>xi + 4>i) + v^i^^J^COS (4 >x2 + ^ 2 ) \ 

cos ( 0^1 + 4 > 1 ) + y/ 527 x 2 COS ( 0 x 2 + ^ 2 ) 

\J ^10yl COS (4>yi + 4>i) + 520y2COs(4>y2 + ^>2) 

\ V^^lTyl cos (0yi + 4>i) + Y/527y2COS (0y2 + ^>2) / 

( 2 ) 

Where 61,62 are constants for a particle. Px 2 ^ 

7ccl5 7 cc25 f^y2'i 0ccl5 4^x2i ^x2’) 0^1? 0^/25 ^yl") ^2/2 

are Ripken’s lattice functions and phase advances. 0i, 02 
are initial phases of a particle. These lattice functions 
at different position of ring have different values, and 
they are periodic functions with the period equal to the 
circumference of the ring, just like the Courant-Snyder 
parameters in uncoupled case. 

1 > Ripken’s lattice functions for CSRm 

Ripken’s lattice functions at different position can 
be calculated from the eigenvectors of the whole ring 
transmission matrixes[5]; Figure 7, 8 show the Ripken’s 
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P functions Px 2 ^ Pyi^ Py 2 , of the ring CSRm and the 
Courant-Snyder parameters /Sy. 



Fig. 7. Ripken’s (3 functions with coupling and the 
Courant-Snyder /3 function without coupling in x 
direction for CSRm. 



Fig. 8. Ripken’s /3 functions with coupling and the 
Courant-Snyder {3 function without coupling in y 
direction for CSRm. 

It can be found that(3(3x2 have the same shape 
with Courant-Snyder parameters/3a,, and Pyi^l3y2 have 
the same shape with/3^. 

2 > The boundary function of coupled beam’s projection 

After obtained the Ripken’s lattice functions of 
CSRm, now the Equation (2) can be used to explain 
beam’s coupled betatron motion in CSRm. From Equa¬ 
tion (2) it can be found that a particle’s projection in 
(x, x’) phase plane or in (y, y’) phase plane is separately 
determined by two modes, the mode 1 parameterized by 
the Ripken’s lattice functions with subscript 1, and mode 
2 with subscript 2. Each independent mode’s projection 
in phase plane with different phase ^i, or 02 is ellipse. 

Hence a particle’s footprints on phase plane at one 
fixed watching point are determined by the two ellipses, 
as Eigure 9 shows. 



-2.51- ^^^^^ -' 
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Fig. 9. A particle’s projection for 20000 turns in 
one point of the ring(blue point). The yellow and 
green cycles present the projections of two modes. 

In simulation one find that the boundary of a parti¬ 
cle’s projection in a watching point is fixed, now the func¬ 
tion of the boundary ellipse can be obtained from Equa¬ 
tion (2). The boundary ellipse should pass through the 
point and the maximum vertical coordinate 

of the boundary ellipse is equal to X^ax^X' 

are determined by Equation (3), (4), (5). 

^max \/ ^10x1 'x/ ^2 0x2 • (3) 

= \/^l7^^*^xl + \/^27 A^^x2- (4) 

X'max = + \/^27 x2- (5) 

Where axi^(ax 2 are Ripken’s lattice functions, which can 
be calculated by whole ring transmission matrixes [5]. 
Use above conditions, the boundary ellipse’s function can 
be obtained: 

^ = ^maxCOs(t) 

a;'=-X'cos(i) - 

Where t is free parameter. Equations (3),(4), (5), (6) are 
also suit for the (y, y’) phase plane, in that case one only 
need to replace all x in above equations by y. 

Eigure 10 shows a particle’s projection for 2000 turns 
and the corresponding boundary ellipse (The yellow el¬ 
lipse) calculated by Equations (6). 
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Fig. 10. A particle’s footprint (blue point) and its 
boundary ellipse (yellow one) in phase plane. 


Form Figure 10 it can be found that the single par¬ 
ticle’s footprint is surrounded by the boundary ellipse 
perfectly. 

From Equation (2), the boundary ellipse function for 
beam can also be obtained: 




s{t) 


[ x'= -Xb'cos{t)-^Xb'^^J^ 

— Xb'^ sin(t) 


Xbnia,x= ■\/ei/ 3 xl-|-\ 

/^2/^x2 5 

(8) 


/ ^2 

Q ^x 2 5 

Px 2 

(9) 

max ^ £ iTxl S” ^ 

/e27x2- 

(10) 

Where the ?i and ?2 are the values of ei, 62, re¬ 
lated to the beam’s coordinates (x, x’, y, y) by Equation 
(2), with these particular values ?i and ?2 the function 
V^i/^xi + \/^2/^x2 should be maximized. 

Above equations are also suit for the (y, y’) phase 


plane, for y direction one just need to replace all x in 
above equations by y. 

Figure 11, 12 show projection of beam and the bound¬ 
ary ellipses are calculated by Equation (7). Eigure 12 is 
for X direction, and Eigurell is for y direction. 



y [m] 

Fig. 11. Projection of 1000 particles for 2000 turns 
in y direction. The yellow ellipse is the matched 
emittance cycle for initial beam. Red points are 
the initial beam. The green ellipse is the bound¬ 
ary cycle. 
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-0.1 -0.05 0 0.05 0.1 

X [m] 


Fig. 12. Projection of 1000 particles for 2000 turns 
in X direction. The yellow ellipse is the matched 
emittance cycle for initial beam. Red points are 
the initial beam. The green ellipse is the bound¬ 
ary cycle. 

From Figure 11, 12, one can found that the Foot¬ 
prints of the beam in x direction and in y direction are 
surrounded by the boundary ellipses perfectly. It veri¬ 
fies that the functions of the boundary ellipses obtained 
from Equation (2) are correct, and the coupled beam’s 
boundary in phase plane can be expressed by the Riken’s 
lattice function successfully. 

In X direction (Eigure 12), the boundary ellipse is 
very close to the matched emittance cycle of the initial 
beam. But in y direction (Eigure 12), the boundary el¬ 
lipse is much larger than the matched emittance cycle of 
the initial beam. This is due to the initial particles’ hor¬ 
izontal emittance (2007r mm mrad) is much larger than 
the vertical one (SOtt mm mrad). 

Erom now, the coupled motion for a single particle 
and for a beam in CSRm are clearly expressed by Rip¬ 
ken’s parameters. The function of projection’s boundary 
in phase plane (Equation (6), (7)) is also obtained and 
parameterized by Riken’s lattice function. 

Comparing with beam’s projection boundary in 
phase plane of the ring, one may care more about the 
envelops in the real space. So in next section a new 
method to calculate the coupled envelops will be shown. 

4 Coupled beam’s envelope 

Assuming that at the beginning the initial particles 
are matched beam in the uncoupled ring (the solenoids 
are compensated by the compensation solenoids per¬ 
fectly), and then the compensation solenoids are turned 
off suddenly. So after the compensation solenoids are 
turned off, the transverse motion will become coupled. 
The envelope of the beam will be different from the un¬ 
coupled case. 

Eorm Equation (2), it can be found that: for a parti¬ 
cle the maximum value of x is equal to \/si/3^i + £2 
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the maximum value of y is equal to 82 Py 2 at 

a positiou of the riug. 

Heuce to calculate the coupled euvelope for a beam, 
the ?i, ?2 (particular values of Ci, 62 ), with which the 
fuuctiou \/^i/ 3 xi+V^ 2 /^x 2 will reach the maximum, should 
be calculated out at every poiut of the riug. Fortu- 
uately, Ci, 62 are coustauts for a particle wheu particle 
moves aloug riug [5], aud the Ripkeu’s lattice fuuctious 
Pxi Px 2 have the same mouotouicity. So Ci, 62 will have 
same values at differeut positious of the riug to make 
(s)+y/ 52 / 5 x 2 (s) to be maximized. So we just ueed 
to calculate ?i, ?2 at oue poiut , aud theu calculate the 
euvelope of the beam by: 

E^{s) = \Jsil3^i (s) + ? 2 / 3 x 2 (s), ( 11 ) 

Ey (s) = \Jei(3yi (s) + \Je 2 (iy 2 (s). ( 12 ) 

The ?!, £2 are calculated by follows: The values of 
Cl, 62 should satisfy equation (13), and at the same time 
they should let function V^i, 5 xi+-\/e 2 / 3 x 2 reach maximum. 

\/ei/3xi cos ($xi + ^ 1 ) + \/e 2 / 3 x 2 cos ($„2 + $ 2 ) = x 

x/ei 7 xi cos ( 0 x 1 + ^ 1 ) + x/e 27 x 2 cos ( 0^2 + ^ 2 ) = x' 

^ y£)^cos($yi + $i) + ye^^cos($y 2 + $ 2 )=y 
cos ( 0 yi + $ 1 ) + cos ( 0 y 2 + $ 2 ) = y' 

7xX^ + 2 q;xXx' + /3xx'^ < 

, 7 yy^ + 2 Q;yyy' + < Ay 

(13) 

Where the 7 ^,, 7 ^, OLy^ /3x, Py are uucoupled twiss 
parameters of a riug. 

By uumerical calculatiou, the ?i, ?2 are ob- 

taiued for CSRm: ?i = 0.00063964077389 aud ?2 = 
0.00009085803501. 
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Fig. 13. The coupled envelope in x direction. 

Using equation (11), ( 12 ) and the Ripken’s lattice 
functions Px 2 , Py 2 of the ring CSRm, the cou¬ 
pled envelope of the beam can be calculated out. Figure 


(13), (14) show the envelopes in x direction and y direc¬ 
tion. To check the reliability of our method, the results 
obtained by beam tracking are also shown in Figure 13, 
14. 



Fig. 14. The coupled envelope in y direction. 

The coupled envelopes of the CSRm calculated by 
the method mentioned above are presented in Figurel3, 
14 by the red lines, and the envelopes in some positions 
of the ring obtained by tracking are shown by the cycle 
marker. It can be found that the tracking results of the 
envelope are accordance with the results calculated by 
numerical method, it verifies that the new method used 
to calculated the coupled beam envelopes is reliable. 

As a compare, uncoupled envelope of the beam also 
be presented by blue line. From Figure 13, it can be 
found that in x direction the uncoupled and coupled en¬ 
velopes are nearly have the same value, the blue line are 
nearly covered by the red line. From Figure 14, one can 
find that in y direction the coupled envelope is much 
larger than the uncoupled one. 

The coupled beam envelope is exceeding the dimen¬ 
sion of aperture (±30mm in dipole, ±50mm in quadru¬ 
ple )[1] in y direction. So if there is no compensation 
solenoid, the coupling induced by cooler’s solenoids will 
cause significant beam loss, and most particles lost in y 
direction. 

5 Conclusions 

The coupled transverse motion induced by cooler’s 
solenoids of CSRm is calculated. To explain the simula¬ 
tion results the Ripken’s approach that used to param¬ 
eterize the coupled betatron motion is introduced. The 
Ripken’s lattice functions are calculated for CSRm. The 
functions of boundary ellipse in phase plane for coupled 
beam are obtained. The coupled envelopes of CSRm are 
calculated by a new method. The conclusions are: 1. If 
there is no compensation solenoid, the coupling induced 
by cooler’s solenoids will cause significant beam loss, and 
most particles lost in y direction. 2. The new method 
used to calculated the coupled envelopes are available. 
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